Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





m 
O 
O 

(N 
Sh 

Oh 
< 

(N 
> 
OS 

o 



(N 
O 



Systematic effects in the measurement of polarization by the 

PLANCK telescope 

G. Franco 1 , P. Fosalba 1 ' 2 , J. A. Tauber 1 

1 Research and Science Support Department of ESA, ESTEC, P.O. Box 299, NL-2200 AG Noordwijk, The 
Netherlands 

2 Institut d'Astrophysique de Paris, 98bis Boulevard Arago, F-75014 Paris, France 
Received 2 October 2002 / Accepted 6 March 2003 

Abstract. We present estimates of the response to polarized signals by the PLANCK telescope. These estimates are 
based on a set of simulations, using a physical optics code (GRASP8) , for linearly polarized detectors at different 
frequencies and located in different positions of an early design of the PLANCK focal plane. We show how the 
optical aberration introduced by the studied focal plane configuration affects absolute and relative orientation of 
the sky signals polarization planes. In addition, we compute the spurious signal introduced by the telescope optics 
into a PLANCK-type measurement of the sky polarization. Our analysis shows that the spurious polarization 
expected in a PLANCK-like experiment is typically of the order of 0.2% of the unpolarized intensity. 
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1. Introduction 

4 : 

The CMB temperature anisotropy carries an essential information on the origin and evolution of large scale structures 
Q. ■ in the universe. Its accurate measurement by experiments probing large (see e.g., IBennett et al. 1994J1 and small 
| scales (see e.g.. lScott et al. 2002l and references therein) already allows us to constrain cosmological models with good 
^ precision (see e.g.. ILewis &: Bridle 20 02). In addition, the quadrupole temperature anisotropy at the last scattering 
. . | surface generates a polarization signal (of order 10% of the CMB intensity) carried by the scalar and tensor modes 
^ ■ of the CMB anisotropy IjRees 1968|l . This polarized anisotropy is a unique probe of primordial gravitational waves 
\ ( |Crittenden et al. 1993a, b| ) and the reionization era of the universe (see e.g., |Zaldarriaga 1997| IHu 2000(1 . Moreover, 
$_j ■ CMB polarization helps break the degeneracy affecting the cosmological parameter estimation from CMB anisotropy 
5^ \ data alone (e.g., |Zaldarriaga, Spergel, ik. Seljak 1997| ). 

Detecting CMB polarization is a serious experimental challenge which requires an unprecedented experimental 
sensitivity, good knowledge and clean removal of foreground contamination, and a careful analysis of systematic 
effects, the latter recently attracting increasing attention (see e.g.. ICarretti et al. 20011 |Kaplan fc Delabrouille 2001| 
|Leahy et al. 200l jj^£gk enko 20011 iFosalba et al. 2002|). Current upper limits for CMB polarizations are at the level 
of 10 uK IjHedman et al. 20011 [Keating et al. 2001] Ide Oliveira-Costa et al. 2002J) . Despite great experimental efforts, 
CMB polarization has remained undetected until very recently with the first detection of -E-mode polarization by the 
DASI experiment tK ovac et al. 2002|l . A new generation of CMB experiments is being specifically designed to image 
the polarization anisotropies in the microwave sky with ground-based telescopes (e.g., COMPASS), interferometers 
(e.g., CBI), balloons (MAXIPOL, B2K), and satellites (MAP, PLANCK). In particular, the PLANCK experiment 1 
will provide multi-frequency full-sky maps of the polarized sky with a sensitivity better than 5/iK for practically all 
the channels IjTauber 20 00 ) . Ultimately, the polarization performance of PLANCK depends on the ability to minimize 
spurious instrumental effects which contaminate the sky signal. 

In this paper, we investigate the impact of one specific instrumental systematic - the optical aberration in the 
PLANCK telescope - in the measurement of polarized signals. The PLANCK telescope is an off-axis aplanatic system 
with a very wide field of view of ~ 8° on the sky. The focal plane is filled with individual corrugated horns that 



PLANCK is the third medium-sized mission in the Horizons 2000 Scientific Program of the European Space Agency, 
http : / / astro . est ec . esa.nl/Planck| 
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collect radiation from the telescope and deliver it to band-limited detectors (with bandwidth ~ 25%) within the 30- 
850 GHz frequency range. Detectors collect largely linearly polarized radiation, defined either by a waveguide hybrid 
discriminator or by a linear grid. This experimental scheme is common to other current CMB polarization experiments 
(e.g., MAP). 

Although PLANCK was not originally devised as a CMB polarization imager, its current design incorporates 7 
polarized channels between 30 and 353 GHz, in the Low- and High-frequency instruments (LFI and HFI hereafter), 
including a total of 88 linearly polarized detectors. PLANCK, to be launched in 2007, is currently in its Phase B 
(detailed design) during which many features, such as the focal plane layout, are being optimised for the polarization 
measurement. The scheme exploits the fact that many detectors in the focal plane are sensitive to the same frequency 
band and are spatially arranged such that their polarization planes are rotated with respect to each other by multiples 
of 45°. In this way, by combining samples from several linearly polarized detectors, it is possible to extract all the 
polarization information in terms of the Stokes parameters (see e.g-. IKraus 19821 ICouchot et al. 1 999). However, an 
optimal measurement relies critically on the similarity of the polarized angular responses, also called beam or radiation 
patterns, of the differently oriented polarized horns. These radiation patterns will be difficult to measure accurately 
on the ground and they will have to be reconstructed in flight using celestial sources such as planets and the Galactic 
plane. 

In this paper we shall study two basic issues regarding the ability of PLANCK to measure polarization: 

— What is the effect of the feed and telescope optics on the properties of the measured polarized sky signal? In 
particular, are the absolute and relative orientation of polarization planes in the sky preserved by the optics? 

— What is the impact of the instrumental optics on the differential measurement of polarization? 

In what follows we shall use Physical Optics (PO) modeling to simulate the interaction of radiation with PLANCK 
telescope optical system. Unlike Geometrical optics (GO), which is characterized by the neglect of the wavelength, 
PO uses the concepts of electromagnetic transverse waves and their propagation properties to simulate the behaviour 
of light in optical systems, being the most accurate method to simulate diffraction, interference and polarization 
phenomena. In this study we use GRASP8, a software package developed by TICRA (1997)[ to compute the polarized 



radiation pattern on the sky for different transmitting feeds located in the focal plane^ (see ijA.ll for further details). 
In particular, we shall simulate the feeds located at three representative (off-axis) positions and at two different 
frequencies, and examine the effect of the optical asymmetries on the measurement of the Stokes parameters. Using 
such simulations we aim at answering the two basic questions mentioned above and draw some general conclusions 
regarding the impact of the PLANCK design under study in this paper, for the measurement of CMB polarization. 

The paper is organized as follows: in we address the problem of simulating and comparing the beam (or 
radiation) patterns generated using PO modeling. Results of the simulations are presented and discussed in and 
our final conclusions are given in In the Appendix we give a full description of the simulations procedures with 
GRASP8 code. 

2. Modeling Far-field Radiation Patterns 

2.1. Introducing reference frames 

Polarized radiation can be fully described by the projection of the electric field vector, E, in a 2-dimensional orthogonal 
basis. For the far field radiation pattern, Ef ar , we shall adopt Ludwig's 3rd definition of co- and cross-polarization 
ULudwig 19731 , 

Efar = EcoUco + E 

cr (7cr (i) 

where a co and a cr are the co- and cross-polarization unit vectors which relate to the usual spherical polar basis (erg, cr^) 
as follows, 

G co = eg sin (p + er<f, cos <j> 

o'er = o~ 9 cos (f) — (j0 sin (f). (2) 
By projecting onto the Cartesian coordinate system, 
erg = er x cos 9 cos cj> + er y cos 9 sin ej> — er z sin 9 

<T0 = — a x sin0 + cry cose/) (3) 

2 By the principle of reciprocity, modeling a feed transmitting from the focal plane to the sky is equivalent to modeling a 
detector horn, in the same position, collecting radiation from the sky, since, by time reversal, the transmission pattern is the 
same as the reception pattern. 
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one finds that at boresight, i.e., the north pole of the sphere (8 — » 0°), it holds that a co — > a y and a cr — > ct^ (see eqs. 

and J2J), so the chosen polarization basis lies along the Cartesian one. Away from boresight, Ludwig's polarization 
basis is obtained by parallel-transporting the local Cartesian basis (cr x , o~ y ) along great circles over the sphere. 



2.2. Defining cross-polarization for non-ideal antennas 

Here we shall define a robust way of characterizing cross-polarization in non-ideal optical systems. In particular, we 
discuss cross-polarization definitions for the case of feeds with arbitrary orientations in the focal plane. 

Let us denote the orientation of a given feed (with respect to its symmetry or rotation axis) by ipFP- We define 
a reference orientation, denoted by ip F r P = 0°, with respect to which we can define any arbitrary orientation of a 
given feed. Similarly, we must introduce a corresponding far-field reference frame with respect to which the sky beam 
patterns from a given feed with an arbitrary orientation in the focal plane, will be referred to. 

For an ideal telescope 3 with an on-axis feed, one can always find a reference frame in the far field for which the 
cross-polarization power (see eq.QJ) vanishes (E^ r — 0). This constitutes a convenient reference frame for sky beam 
patterns of a given feed. However, for a real telescope, any asymmetry in the antenna system will generally introduce 
a non-vanishing cross-polarization component in the far-field radiation pattern. Moreover, non-ideal telescope optics 
also introduce a differential rotation of the principal plane of polarization in the far field with respect to that physically 
applied to the feed in the focal plane. 

We recall that a proper measurement of polarization in terms of the Stokes parameters of the sky signal requires 
that the orientation of the feeds in the focal plane be such that the principal planes of polarization on the sky are 
spaced by 45°. Therefore, it is important to choose an appropriate method for determining this principal plane of 
polarization, which is here defined as that in which the Cross Polarization Discrimination (XPD) is maximized. 

We define maximal cross-polar discrimination angle, (/>xpd, as the angle by which one should rotate the reference 
frame in the sky in order to maximize the Cross Polar Discrimination, XPD. Below we introduce three different 
methods for the computation of 4>xpd, each one relying on a different definition of XPD 4 : 

— The 1st Method is the most commonly used in the literature. 4>xpd is the angle that maximizes XPD defined 
as the ratio of the (single) co-polar power peak to the highest cross-polarization power peak. In dB units we can 
write the maximized XPD as, 



AW max = AW(^xpd) = 10 lo gl 



Max[ElfEl 



Max{E r c ° t E r c ° t *) 



10 log 10 [Max( {Elff )] - 10 log 10 [Max{ [E^f )] (4) 
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In the 2nd Method, XPD is defined as the difference between the co-polar power peak and the cross-polarization 
power along the same direction. Since (frxPD is the angle that maximizes this difference, we can write 

AW max = AW^xpd) = 10 log 10 [Max( (£™*) 2 )] - 10 log 10 [(E r c f{6 Max , <t>Ma*) ?] (5) 

where 8 Max and 4>Max are the coordinates of the co-polar power peak in the spherical coordinate system [8, (f>). 
The 3rd Method determines <Pxpd as the angle that maximizes the difference between the integrated co-polar 
and cross-polarization amplitudes over the main beam, and thus the maximized XPD is given by 

(6) 

where the integration interval on 8, [—8q; +9q], must be taken wide enough to encompass most of the power in the 
co- and cross-polarization patterns. We chose 8q to be at least twice the Full Width at Half Maximum (FWHM) 
of the co-polar power pattern. 

Note that in the equations above E™* and E™* are the (complex) co- and cross-polarization components of the 
radiation pattern, after having performed a rotation of the far field coordinate system by an angle (4>xpd in this case) 
with respect to the original reference frame (i.e, the far-field frame defined for a non-rotated feed in the focal plane), 

E r c ° l = E cr sin cj> X PD + E co cos <j> X PD 

Elf - E cr cos (j) X PD - E co sin (j) X PD (7) 

In general, the far-field radiation pattern of an asymmetric radiating system fed by a corrugated horn has a 
cross-polarization pattern with a multiple peak (rather than a simple peak) structure with several maxima which are 

3 We define an ideal telescope as one with rotationally symmetric reflectors and on-axis fully linearly polarized detectors. 

4 For simplicity of notation, XPD will always be expressed in dB units as a power difference, AW. 
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not spatially coincident with nor symmetrically placed with respect to the co-polar peak (e.g. see Figure lA~4|l . Such 
an asymmetric cross-polarization pattern for non-ideal antennas makes the cross polar discrimination angle defined 
in the 1st method potentially ambiguous. In principle, the 3rd method should be appropriate for measurements of 
diffuse polarized emission, as it is based on integrating radiation filling the main beam (i.e, an extended region around 
boresight), while the 2nd method is more appropriate for measuring polarization from point sources as it relies on a 
local measurement of power radiated. 



2.3. Differential rotation of the polarization plane 

Let us consider again an ideal antenna system. In this case, it is expected that a rotation of a given feed around its 
own symmetry axis, from an initial position, ip FP , to a final position, ippp, would result in an equal rotation of the 
principal plane of polarization of the radiation pattern in the sky J . Thus we can write, 

0XPD ~ 0XPD = -bl>FP-^Fp) ( 8 ) 

Setting <Pxpd = fi'xPD as the maximal cross-polarization discrimination angle for the beam pattern when the feed is 
in its original or non-rotated position (V->pp = 0°), the previous equation becomes, 

4>XPD ~ <$>XPD = ~^FP ■ (9) 

This equation holds for any ideal antenna system if the polarization properties of the radiation pattern are equally 
preserved for all orientations of the feed in the focal plane. 

However, for a non-ideal antenna system, such isotropy is broken and, due to factors such as the asymmetry of 
the reflectors or the off-axis positioning of the feed, equation will only hold to a first-order approximation, i.e., 
4>xpd — ftxPD ~ ~i>FP- Consequently, in order to quantify how good this approximation is (or on the contrary, how 
non- ideal is the antenna), we shall define the residual angle, A^, as the additional angle to which the reference frame 
of the beam pattern in the far-field should be rotated, in addition to (—ippp), so as to have the co-polar axis of the 
beam pattern reference frame aligned with the principal polarization direction. Thus we define 

A = (j>XPD - 4>X PD - (-TpFp) (10) 

where 4>xpd is computed using the methods described in fl2.2l and 4> r xpD ^ s the angle <Pxpd for the non-rotated feed, 

4>XPD =<t>XPD(lpFP = 0°). 



3. Simulation Results 

In what follows, we shall focus on determining this residual angle, along with the associated angle 4>xpd and power 
difference /SW(4>xpd), making use of the PO simulations described in flA.ll In the present study we compute the sky 
beam patterns for PLANCK-LFI 30 and 100 GHz feeds at different positions in the focal plane unit, for a range of 
different orientations of the feed. The orientations studied correspond to rotations around the symmetry axis of the 
feed by angles 6 ipFP = 0°, 45°, 60°, 90°, 120°, 135° and 180°. In particular, we carry out numerical simulations of the 
far-field patterns fed by corrugated horns located at three different positions in the focal plane ranging from very close 
to the optical axis, to the edge of the field of view (see Figure lA~3|) . 

This will allow us to draw some general conclusions about the effect of the feed positioning on the preservation of the 
polarization properties through the antenna. We emphasize that, despite concentrating on the PLANCK experiment, 
our results can be useful for other mm/cm-wave polarimetric experiments with similar optical systems. 



3.1. Effect of feed position and orientation 

Tables n 121 an( i El display the values of the residual angle (A^) for all the simulated cases. These values were computed 
using equation (|1UI) . The corresponding values for 4>xpd and AW max were determined using the three different methods 
described in i)2.2l In each table, each A^-column shows the results for a different location of the feed in the focal plane 
(see Figure lA~3ll . We point out that, when using the 3rd Method, the main beam pattern is integrated over a solid 
angle centered at boresight (the north pole of the sphere) and down to a latitude 9q = 1.0° and 0.5° for the 30 and 
100 GHz feeds, respectively. These solid angles are chosen so as to include, at least, all radiated power down to the 
FWHM (see ggH) . 

5 In fact, in a dual reflector system such as that of PLANCK, the rotation of the principal plane of polarization on the sky 
would be symmetric and not equal to the rotation of the feed, due to the "mirroring" of the radiation pattern in the beam 
coordinate system when the signal is reflected from one mirror to the other. Equation (JHJ takes this symmetry into account. 

6 Note that, by symmetry, a feed rotation by ipFP is equivalent to a ipFP — 180° rotation, so we need not consider cases with 

Tp FP > 180°. 
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Results on the Determination of A^ by Maximizing the Difference 
Between the Co- and Cross-polar Power Peaks (1st Method) 



30 GHz 




Position 1 




Position 4 


Position 27 


4>FP 


<j>XPD 


AW Max 




4>XPD 


AW Max 


A^ 


4>XPD 


AW M ax 
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n 
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(dB) 


n 


n 


(dB) 


n 
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(dB) 


n 





0.4 


33 





4.8 


29 





7.9 


26 





45 


-43.9 


33 


0.7 


-39.1 


29 


1.2 


-35.1 


26 


2 


60 
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33 


0.7 
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29 


1.2 
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26 


2.4 


90 


-89.3 
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0.3 
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29 
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26 
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29 


1.0 


-110.8 


26 


1.3 
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-135.0 


33 


-0.4 


-129.8 


29 


0.5 


-126.5 


26 


0.6 


180 


-179.6 


33 
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29 





-172.1 
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100 GHz 
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Position 4 
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AW Max 


A 
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0.3 


32 
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30 





9.3 


28 


0.0 
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0.9 
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-114.8 
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-110.4 
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135 
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32 
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-129.5 


30 
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-125.6 


27 


0.1 


180 


-179.7 
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-174.1 


30 





-170.7 


28 






Table 1. Maximal cross-polarization discrimination angle 4>xpd, AWmax and the residual angle A^ for the 30 GHz and the 
100 GHz feeds on positions 1, 4 and 27, for different orientations tpFP of the feeds being cfixPD determined by the 1st Method 
described in H2.2l and AWMax evaluated using equation I^J. 



Results on the Determination of A^ by Maximizing the Difference Between 
the Co-polar Power Peak and the Cross-polar Component (2nd Method) 



30 GHz 
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100 GHz 
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Table 2. Same as Tabled with 4>xpd being determined by the 2nd Method, as described in H2.2I 



Our results show that the estimated differential rotation of the polarization plane (see H2.3|) in terms of the 
residual angle, A^, depends strongly on the definition adopted for maximum cross-polarization discrimination and the 
associated angle, 4>xpd- In practice, we see from the calculations that the 3rd Method gives what one would intuitively 
expect from a robust estimate of cross-polarization: there is no differential rotation of the polarization plane in any 
of the cases studied (i.e., the residual angle, is found to be compatible with zero in all cases, see TableOJ). In fact, 
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Results on the determination of by Maximizing the Difference 
Between the Integrated Co- and Cross-polar Power Beams (3rd Method) 



30 GHz 
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-54.4 


28 





-51.3 


25 





90 


-89.5 


32 





-84.4 


28 





-81.3 


25 





120 


-119.5 


32 





-114.4 


28 





-111.3 


25 





135 


-134.5 


32 





-129.4 


28 





-126.3 


25 





180 


-179.5 


32 





-174.4 


28 





-171.3 


25 







100 GHz 
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45 
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33 





-39.4 


29 


0.1 


-36.2 


26 


0.1 


60 


-59.5 


33 





-54.4 


29 


0.1 


-51.2 


26 


0.1 


90 


-89.5 


33 





-84.4 


29 


0.1 


-81.2 


26 


0.1 


120 


-119.5 


33 





-114.5 


29 





-111.3 


26 





135 


-134.5 


33 





-129.5 


29 





-126.3 


26 





180 
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33 





-174.5 


29 





-171.3 
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Table 3. Same as Tabled with 4>xpd being determined by the 3rd Method, as described in 32.21 



the maximum value of the differential rotation found (A^ = 0.1°) for the 100 GHz feed at the off-axis positions 4 
and 27 (see Figure |A.3|> . is not significant given the errors, typically of this order, coming from the accuracy in the 
computation of the beam pattern itself as well as its finite sampling in spherical cuts. On the other hand, Methods 1 
and 2 are misleading when measuring diffuse polarization and they indicate a significant rotation of the polarization 
plane, whose magnitude varies systematically with frequency and position in the focal plane (see Tables ^ an( l El ■ 

In summary, we see that the 3rd Method turns out to be an appropriate definition of cross-polarization and the 
related cross-polar rejection angle for measurements of diffuse polarized emission. Furthermore, making use of this 
definition, we conclude that the PLANCK telescope does not introduce any second order (i.e. detector orientation 
dependent) spurious polarization. 

For the absolute orientation of the polarization plane of the sky signal there is a misalignment angle introduced 
by the telescope optics. Referring to Figures TA . 1 1 and f A . 21 by misalignment angle we mean the angle, as seen from the 
the xy-plane of the line-of-sight Cartesian coordinate system of the telescope 7 (Cpm), between the y-axis of the feed 
coordinate system (e.g. fl) and the y-axis of the corresponding coordinate system of the sky beam pattern (mbl as 
corresponding to fl). 



Feed Position in Focal Plane 


1 


4 


27 


Misalignment Angle 


0.3° 


1.9° 


4.3° 



Table 4. Misalignment angle between the y-axis of both the feed and the far-field coordinate systems, for each position of the 
feed in the focal plane. Feed positions are shown in Figure lA~3l 



Table 0| shows the results for the misalignment angles, computed for each simulated position of the feed in the focal 
plane. They take into account the values of maximal cross-polarization rejection for the non-rotated feed, (ipjrp = 0°) 
(as determined using the 3rd Method). It can be seen that the bias or misalignment angle of the polarization plane 
increases as the feed moves from the center (on-axis) towards the edge of the focal plane (as we would expect), yielding 
an upper limit of about 5° . Note that these results are independent of feed orientation and frequency, within the errors 
(see Table 0J, which means that they do not affect the differential measurement of polarization for each position of the 
feed. In addition, this bias can always be corrected by rearranging the orientation of each detector in the focal plane. 

7 By definition, the z-axis of the line-of sight coordinate system points to the direction of maximum power radiated (or 
received) by an on-axis detector. 
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3.2. Predicting far- field beam patterns 

PO modeling of a large telescope such as PLANCK requires computer intensive calculations with processing time 
scaling as the fourth power of frequency. Moreover, modeling the response of a focal plane array including many 
detectors at several frequencies becomes an extremely demanding task. Therefore, we investigate to what extent it is 
possible to use a known (pre-computed) model of the main beam response (i.e, the peak of detector angular response), 
for a given feed and focal plane layout, to estimate or predict the sky main beam patterns for that detector in 
different focal plane configurations, i.e., for the feed in the same off-axis position but different polarization sensitivity 
directions. In order to do so, we must assess the similarity between the main beam patterns for different orientations 
of the radiating feed in the focal plane. At each frequency, we compare the beam pattern generated from two different 
feed orientations for which the far- field reference frame has been rotated by 4>xpd (with its co-polar component aligned 
with the principal plane of polarization). Figure ^ shows the differences (in dB units) between the power contours of 
several pairs of main beam patterns. For illustrative purposes, we display the difference contours for the case where 
the definitions of maximum cross-polarization discrimination discussed in A2.2I disagree the most (Figures ^a)^b)) 
and for the cases where the feeds are closest to the center of the focal plane (Figures IHc)QId)). 

As can be seen from Figure^ co-polar power differences are always below 3.0 dB, except for a few point-like regions, 
while cross-polar power mismatches are locally found up to 15 dB. Since cross-polarization peaks are typically about 
30 dB bellow the co-polar peak (see AWmux values in TableOJ, it is sufficient to concentrate on the dominant (co-polar) 
component of the polarized beam pattern to estimate how much the beam pattern changes as a function of the feed 
orientation. No significant difference is found between the co-polar patterns on Figures Gfa) and^b) and therefore we 
conclude that these difference contours cannot be used to determine a robust definition for cross-polarization. 

More interestingly, from Figure ^ we conclude that using a pre-computed sky main beam pattern from a feed at a 
given frequency and orientation in the focal plane (and the corresponding maximal cross-polarization rejection angle, 
4>xpd), it is possible to predict the shape of the main beam pattern in the sky for any other orientation of the same 
feed, within ^3 dB (~15 dB) for the the co-polar (cross-polar) pattern. 



3.3. Spurious optical polarization 

The presence of cross-polarization and the differences in the shape of the polarized beams patterns from a given feed 
with different orientations around its symmetry axis will result in an observation of a polarized signal not present in 
the sky, but rather generated by the instrument asymmetries. This is what we shall call "spurious optical polarization" . 
We quantify this effect in the modeled telescope as the response of the system to a diffuse unpolarized sky signal, in 
terms of normalized Stokes parameters. Following standard conventions l|Kraus 1982(1 . PLANCK will measure Stokes 
parameters I, Q and U by combining the output of four linearly polarized detectors; a first pair having principal planes 
of polarization at right angles to each other and a second pair with principal planes of polarization at right angles as 
well, but rotated by 45° with respect to the first pair of detectors. For this purpose we shall use the simulations results 
for the 30 and 100 GHz feeds at different positions in the focal plane unit and orientations 4>fp = 0°,45°,90°, 135°. 
The power available to the detectors is proportional to the convolution of the beam with the sky signal, i.e. 



Poc f dn(s B -s s ) 



(11) 



where Sb = (Ib,Qb, Ub,Vb) is the Stokes vector of the beam on the sky and Ss is that of the sky signal. If the sky 
signal is uniform, then 

Poc / dn (E 2 co + El) (12) 

We shall use normalized Stokes parameters, qs = Qb/Ib, and ub = Ub/Ib, to characterize the spurious optical 
polarization, 

_ P(^) - P(y>! + f ) _ P(V; 2 )-P(^ + f) 

qB P^O+P^i + f) ' UB P^x)+P{^ + 1) { ' 

where P(ipi) + P^Pi + § ) = P{4 } 2) + P(^2 + f ) an d P(V0 is the total power integrated over the beam (within a 
co-latitude range [— 8 ; +9q\). V^2 — "01 = 45° is the relative phase shift in the polarization planes. The calculated 
estimates of those parameters for the modeled feeds are shown in Table for different limits of integration da- It can 
be seen that qB and ub depend on the chosen integration limit, if it less than a few times the FWHM. 

Table [S] also shows that the spurious polarization depends strongly on location and frequency and that a high 
level of up to 5% polarized contamination of the diffuse unpolarized intensity is found for feeds at 30 and 100 GHz. 
The observed large variation of this level across the focal plane seems, at first glance, difficult to explain. In order to 
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0.5° 3.0° 
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-0.275 
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3.17 
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0.00117 


0.00110 


5.2 


5.1 


5.0 


-0.427 -0.204 


-0.205 




100 GHz 




Position 1 




Position 4 


Position 27 


00 


0.2° 


3.0° 


20.0° 


0.2° 


3.0° 


20.0° 


0.2° 3.0° 


20.0° 


IB (%) 


0.164 


0.093 


0.093 


-0.138 


-0.0398 


-0.0422 


1.42 2.05 


2.05 


U B (%) 


-0.00418 


0.0052 


0.0052 


-0.51 


-0.086 


-0.086 


2.22 3.44 


3.43 



Table 5. Normalized Stokes parameters gs and ub for different limits of integration 9o; at FWHM, 3° and 20°. These 
parameters represent the optical spurious polarization introduced by PLANCK telescope optics for the studied focal plane 
configuration (in bold) and other study cases (see text). 



illustrate the issue, we show in Figure |21 the dependence of qs and ub across the main beam for the two worst cases 
in Table El the 30 GHz feed on position 4 (Fig. [^a)) and the 100 GHz feed on position 27 (Fig. Hfc)) and for the 
cases where the two feeds are in their actual positions in the studied PLANCK focal plane layout; the 30 GHz feed 
on position 27 (Fig.[3») and the 100 GHz feed on position 4 (Fig.^d)). 

It can be seen that qs and ub exhibit large fluctuations across the main beam, which may cancel out more 
effectively when integrated in angular space. The final integrated values of qs and ub will therefore depend on the 
fine details of the modeled patterns, which are affected by the exact location of the feed in the focal plane. We can 
expect that, in a real system where the detectors are sensitive to wide bandwidths, part of the fluctuations will be 
averaged out. In addition, it is not clear to what extent the fine details of the models will be reproduced by the real 
physical systems. For these reasons we can consider that the levels of spurious optical polarization here estimated are 
very conservative upper limits. 

It is however emphasized that, given the expected small relative amplitude of the polarized with respect to the 
unpolarized CMB signal (about 10%), this upper limit of optically-induced polarization cannot be ignored. At these 
frequencies, there are two large scale sources of unpolarized diffuse components in the sky: the isotropic CMB emission 
(at T ~ 2.73 K) and the CMB dipole, a Doppler temperature anisotropy in the sky of amplitude ~ 3mK (see, e.g. 
ILineweaver et al~9 96 ) . As PLANCK scans the sky, spurious optical polarization will introduce a constant bias for 
the former and a slowly varying bias for the latter. The level of these biases (~ 30toK and 0.03mK, for each 1% of 
spurious polarization) is quite large compared to the expected CMB polarized signal of ~ 5 — 10^/K. Note however 
that a constant bias does not affect the measurements of small scale anisotropics, and a slowly varying bias is likely 
to be easily removable using calibration procedures. 

The main source of concern in the detection of polarized sky signals for a PLANCK-like experiment will arise from 
the spurious signal that leaks from the measured unpolarized intensity itself. The percentage of spurious polarization 
is, in the very worst cases, only factors of a few smaller than the percentage of polarization signal. Potentially this 
spurious level will ultimately introduce a fundamental limit on the measurement of CMB polarization anisotropics at 
sub-degree scales, i.e. at angular scales comparable to the main-beam size. It is however emphasized that the spurious 
polarization estimated for the actual positions of the PLANCK feeds is much lower than in the worst case scenario, 
and it is certainly below the expected polarized CMB signal. Moreover, this situation should be further improved if we 
take into account other relevant effects such as the large detector bandwidth that tends to suppress to a large extent 
the small-scale angular variations observed on monochromatic radiation patterns. 

3.4. Model uncertainties 

This study has not taken into account imperfections in the optics due to the effects of dust contamination, microcracks, 
asymmetry of the focal plane, or the effects of optical misalignment, since it is today very difficult to model such effects. 
However, qualitatively it can be argued that: 

— As long as dust grains are not elongated and preferentially aligned, their effect is largely to scatter power from the 
main beam into the far sidelobes. In this study we are interested in the polarization introduced through the main 
beam and therefore dust can be safely ignored. 

— The effect of microcracks is qualitatively similar to that of dust with the exception that elongation and preferential 
alignment are expected. Fortunately, in the case of PLANCK, all technological indications are that no microcracks 
will be present. 

— Misalignment introduces additional aberration in the optics that translates into distortion of the main beam. 
Distortion is always related to depolarisation and therefore it is of potential concern. However, the aberrations 
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introduced by misalignment are, by design and test, kept at a level which is a small fraction of the intrinsic optical 
distortions. Therefore it is reasonable in this initial study to assume that it will remain a second order effect. 

— An asymmetric focal plane will likely introduce differential effects in the response of a single horn to different 
planes of polarization. However, the effects of mutual coupling between PLANCK horns that are adjacent in the 
focal plane are kept, by design, at a rather low level with respect to the main beam (i.e, at least 30 dB below peak) 
and we assume here that they are not affecting the main results of this study. 

In PLANCK, the higher frequency detector chains include also a number of optical elements (such as filters, detector 
embedding structures, cavities, etc.) which potentially introduce polarization effects. Such chains have been measured 
and found to introduce very low levels of cross-polarization. However, since they are of the same order as the effects 
listed above, they should be included in future simulations. 

Finally, we emphasize that this study is based on a numerical model, which naturally is subject to some degree 
of uncertainty when compared to reality. In this respect, the software used for our simulations (GRASP) is today 
considered to be the benchmark software tool for analysis of large antenna systems. Recent studies | |Murphy et al. 2001 ) 
have been carried out to investigate the reliability of a variety of software tools, including GRASP8, in modelling 
optical systems in the submillimetre wave band. In these studies, simulated results from different software packages 
were compared to each other as well as to experimental measurements. These studies show GRASP to be the most 
reliable package, agreeing with experimental measurements down to a —40 dB level in predicting the amplitude pattern 
and showing the best agreement in phase measurements. Moreover, the same study has shown that GRASP8 is the 
only available package capable of handling polarization in the case of off-axis reflectors systems. 

Therefore, in terms of model uncertainties, it is considered that this study is as close to real as can be expected, 
within the assumptions made. However, it is acknowledged that the real system performance must ultimately be 
measured. A comprehensive pre-launch test campaign is planned for PLANCK; these measurements will be later 
combined with in-flight characterisation to obtain the final system performance. 

4. Summary and Conclusions 

We have studied the systematic effects introduced by the PLANCK optical system in the measurement of polarized 
signals. For the purpose of this study we have carried out PO simulations of the former PLANCK telescope configu- 
ration with the 30 and 100 GHz feeds at three different feed locations (see Figure lA~3)l and several orientations. From 
the analysis of these simulations, our results (see fJ^J) yield the following main conclusions: 

— The relative orientation of polarization planes of diffuse sky signals is near-perfectly preserved by the optics, so 
that no systematic effects in the measurement of polarization are to be expected from the feed orientations in 
the current baseline. This symmetry is only observed when a proper definition of cross-polarization for non-ideal 
antennas is used (see H2.2| l. 

— The absolute orientation of polarization planes is systematically rotated by the telescope. Such a misalignment 
angle can be as large as 5° and depends only on frequency, increasing monotonically as the feed is moved towards 
more off-axis positions (see Table 0J. However, this misalignment can be corrected by appropriate design of the 
focal plane. 

— Due to the high symmetry of the telescope, it is possible to use a pre-computed polarization main-beam pattern 
for a given feed to predict main beam patterns (within 2° in co-latitude) for arbitrary orientations of that feed 
within ^3 dB (~15 dB) for the the co-polar (cross-polar) pattern (see H3.2(l . This will allow systematic studies of 
the polarization measurements, involving a larger number of feeds, by using a minimal amount of CPU time. 

— Spurious optical polarization (i.e, polarization solely introduced by the telescope optics) is estimated to be less than 
0.2% for the actual locations of the 30 and 100 GHz feeds in the studied focal plane of PLANCK (see Figure lA~3)l . 
Other (less optimal) studied cases show spurious optical polarization as large as ^5% of the unpolarized sky signal 
(see Table [5jl . It is noted however that these estimates rely on uncertain fine details in the modeled patterns and 
therefore should be considered as order of magnitude only. It is likely that real experiment features such as broad 
bands, detector responsivity and optical imperfections would affect the results, reducing the magnitude of this 
effect. Hence these simulations should only be used as upper limit guidelines. 

Further studies on the impact of PLANCK optics on the detection of polarization should consider the following 
topics: 

— Use of the final optical configuration of the telescope, in particular the focal plane, has recently undergone some 
modifications. 

— Take into account second order effects (already discussed in t!3.4|) that may be important to the detection of 
polarization in the main beam region. 
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— Extending the studied region further into the sidelobes, to include any second order effects that may have an 
impact on the detection of polarization in this regions. 

Finally the results of these studies will be compared to measurements obtained on qualification models of the 
PLANCK optics during a pre-launch test campaign, planned to start in 2004. 
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Appendix A: GRASP8 Simulations 

A.l. Simulations setup and inputs 

As stated before, the telescope configuration used for the present analysis is an early version of the current PLANCK 
telescope design. It corresponds to a dual reflector optical system with an offset aplanatic geometry. The exact telescope 
configuration is shown in Table lA~Tl 

Figures fA . 1 1 and I A . 21 show two different views of the PLANCK Telescope system, with the 100 GHz feed located 
in three different positions in the focal plane (positions 1, 4 and 27, also shown with more detail in Figure lA~3l) . 

Reference frames fl, f4 and f27 are used to describe the three different positions of the feed in the focal plane 
while reference frames mbl, mb4 and mb27 are used to describe the sky beam patterns, the latter having their z-axis 
(in blue) always pointing to the center of the main beam in the sky, set to be the maximum of the co-polar pattern. 
Coordinate systems fl, f4 and f27 are defined on the basis of the coordinate system located on the vertex of the ellipse 
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Main Reflector 



Surface Shape 
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Off-axis Ellipsoid 
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1641.58mm 
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328.15mm 


Projected Rim Shape 
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Rim Major Axis Projected Aperture 


509.30mm 


Rim Minor Axis Projected Aperture 


390.30mm 



Table A.l. Early PLANCK Telescope configuration used for the present study (current telescope dimensions are slightly 
different). 



that subtends the primary mirror (Vpm) while mbl, mb4 and mb27, are defined on the basis of the coordinate system 
located on the center of the primary mirror (Cpm) which has the same orientation as Vpm (see Figures lA~71 and lA~2|) . 
Furthermore, the orientation of Vpm is the same as the spacecraft coordinate system rotated by +13.75° around its 
y-axis. 

For all simulations in this study, we used the polar coordinates (9, cf>) that can range to cover the whole sky 
(i.e., — 7r < 9 < tt and < <fi < ir) s producing ASCII data files that have a structure such as to hold the values of the 
two complex amplitudes (real and imaginary) of the two components of the beam pattern (co- and cross-polar) , for all 
given sets of points (6, <fi). All simulations were done in order to determine the radiation pattern in the sky area around 
the main beam. To sample this part of the sky, coordinate 9 was chosen to range within the limits —20° < 9 < +20° 
in a set of 401 points with intervals of 0.1° while coordinate cf> ranged within 0° < <p < 180° in a set of 36 polar cuts 
separated by 5°. 

In a GRASP8 simulation, the scattered field can be calculated by using Physical Optics (PO) combined with the 
Physical Theory of Diffraction (PTD) or, alternatively, using Geometrical Optics (GO) combined with the Geometrical 
Theory of Diffraction (GTD). According to |TICRA (1997)| for focused reflector systems PO/PTD is intended to be 



used in the far field around beam maxima, whereas GTD can be used in the side-lobe region. For these reasons and 
since we were only interested in the main beam sky area, our simulations were done with PO/PTD calculations. In 
what concerns the density of the PO integration grid, to determine the number of points for which the PO currents 



are calculated in each reflecting surface, we referred to a formula given by TICRA (1997) where the wavelength, 
the reflector diameter and the maximum latitude, 9 max , are input parameters in the determination of the minimum 
number of points for which the PO integration converges. 

To simulate the conical corrugated horns used in the PLANCK telescope, the input radiation signal of the feeds 
used in our simulations had the form of a spherical wave expansion. Figures lA.51 and IA.4I show the far field radiation 
pattern when the input signal radiates directly from the feed to the sky, without going through the telescope system; 
Figure lA~5l shows the plot of the polar cut at 4> = 20° while Figure lA~4l shows the whole contour in the (U,V) plane. 

Note that to project the polar coordinate system (9, 4>) onto the (U, V) plane we have the following transformation, 

U = sin 9 cos <f> ', < 9 < tt (A 1) 

V = sin 9 sin <j> ; < <ft < tt ^ ' ' 

In Figures FA. 41 and fA. 51 note that the beam maxima, the co-polar peaks, are not centered at ([/, V) = (0, 0) nor at 
latitude 9 — 0°. This off-center positioning of the beam pattern in the coordinate system Cpm when the feeds radiate 
directly to the sky is due to the off-axis positioning of the feeds in this same reference frame. When the radiation 
pattern goes through the telescope system instead of going directly from the feed to the sky, a much smaller off-center 
positioning of the beam can be observed in the coordinate system Cpm- However, in order to have a centered beam 
maximum in all plots and contours, for each of the three different feed positions in the focal plane (1, 4 and 27), the 



Note that 6 is a latitude coordinate, while is a longitude coordinate. 
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coordinate system Cpm was reoriented in order to always have the z-axis pointing to the co-polarization power peak. 
This procedure originated the three different coordinate systems (mbl, mb4 and mb27), already mentioned before 9 . 

A. 2. Simulations outputs and results 

Figure lA~6l shows the far field radiation patterns when the input polarized signal goes through the telescope system. It 
can be seen that, as a result of the far field reference frame reorientation mentioned in the previous section, the co-polar 
peaks are always centered at 9 — 0° whether the feed is on position 27 (left plot) or on position 4 (right plot). We 
show the contour plots of the far field radiation pattern around the main beam for various simulations. In particular, 
Figures lA.71 and IA.8I show radiation patterns in the non-rotated far field reference frame and in the reference frame 
rotated by the 4>xpd angle. This angle is determined by each one of the three different methods described in H2.ll 
Plots are shown side by side for an easier comparison. 

As expected, for simulations with a non-rotated feed, i.e., ipFP — 0° fFigure lA~8j) . the highest peak is at the co-polar 
direction before and after the far field reference frame is rotated by 4>^[p D (being 4> t xpd always a small angle), while 
for simulations with a feed orientation of tpFP = 90° (Figure TA.7JI . the highest peak is at the cross-polar direction 
before the far field reference frame is rotated by 4>xpd which then aligns the highest peak with the co-polar direction. 
On the other hand and also as expected, for simulations with a feed orientation of ipFP — 45°, 135° ("Figure lA~9)l . the 
co- and cross-polar patterns have similar peaks before the rotation of the far field reference frame by <Pxpd- 



9 These far field reference frames (mbl, mb4 and mb27) can be seen in Figures I A . II and I A . 21 Note that they have a common 
origin with Cpm, at the center of the primary mirror, but different orientations. 
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(a) 1st Method 



(c) 3rd Method 



30 GHz Feed on Position 27 
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30 GHz Feed on Position 1 
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(b) 3rd Method 



(d) 3rd Method 
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(y» = o°) - tv- = 90=) 

Rotated Co-polor Components dB 




Rotated Cross— polar Components dB 
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100 GHz Feed on Position 1 
(VVp = 0°) - (i/< F p = 45°) 

Rotated Co-polor Components dB 



-0.0201 

-0.020 -0.010 0.000 0.010 0.020 



Ro:a:ed Cross-polar Components dB 



4» 



-0.020 
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Fig. 1. Differences (in dB units) between main beam patterns in far field reference frames rotated by ipxPD, for feeds at 
different positions and orientations in the focal plane. Cross-polarization definitions and 4>xpd angles are determined 
by Methods 1 & 3 described in 5^2 (see Tables □ and 01 . 
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30 GHz Feed at Position 4 30 GHz Feed at Position 27 

Non-rotated Stokes 11(45/135) % Non-rotated Stokes 11(45/135) % 




Fig. 2. Contour plot of the distribution of the normalized Stokes parameters across the main beam (with |0| < 3°). 
The plot shows one of the parameters (ub or qs) for the two worst study cases in Table El (a) the 30 GHz feed on 
position 4 and (c) the 100 GHz feed on position 27. The same parameters are also shown for the feeds in their real 
positions in the studied PLANCK focal plane layout: (b) the 30 GHz feed on position 27 and (d) the 100 GHz feed 
on position 4. 
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X 




Fig. A.l. Cross-section of the objects that define the Alcatel version of PLANCK telescope in the GRASP8 model 
of our study. Besides the primary (PM) and secondary (SM) mirrors, this plot includes the 100 GHz feed in three 
different positions in the focal plane with the corresponding coordinate systems for the feed positions (fl, f4, f27) and 
the different orientations of the main beam in the sky (mbl, mb4, mb27). It also includes the fixed coordinate systems 
for the vertex (Vpm) and the center (Cpm) of the primary mirror and for the spacecraft (SC). 
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Fig. A. 2. Oblique view of the objects that define the Alcatel version of PLANCK telescope in the GRASP8 input file 
of this study, as shown in the previous Figure IA.1I 




Fig. A. 3. Former design of PLANCK Telescope LFI focal plane showing the detector positions used for this study: 
on the top left-hand corner, positions 27 and 4 with the 30 and 100 GHz detectors, respectively and position 1, closer 
to the center of the focal plane unit, with the HFI 857 GHz detector. 
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30 GHz Feed Radiation Pattern 100 GHz Feed Radiation Pattern 




Fig. A. 4. Contour plots of the far field radiation pattern (\0\ < 30°) when the feeds radiate directly to the sky, i.e., 
the input polarized signal does not go through the telescope system. 



30 GHz Feed Radiation Pattern (<fi = 70°) 100 GHz Feed Radiation Pattern (0 = 20°) 
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Fig. A. 5. Polar cuts, at <j> = 20°, of the contour plots in Figure lA~4l 
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Fig. A. 7. Contour plots of the Co- and Cross-polar patterns around the main beam (\8\ < 2°) for the 30 GHz feed 
rotated by ipFP = 90° at position 27. The figure shows (a) the beam pattern in the non-rotated far field reference 
frame and the same beam pattern in the far field reference frame rotated by (b) 4>xpd = —79.6° (see table QJ, (c) 
4 > xpd = -81.5° (see table^ and (d) cj> X PD = -81.3° (see table EJ. 
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100GHz Feed on Position 4 
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Fig. A. 8. Same as Figure \Kj\ (with \6\ < 1°) for the 100 GHz non-rotated feed {ip FP = 0°) at position 4. In this 



figure the far field reference frame is rotated by (b) 4>xpd — 5.9° (see tabled, (c) 
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5.5° (see tabic EJ. 
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6.1° (see table and (d) 
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Fig. A. 9. Same as Figures and El for the 100 GHz feed rotated by 4> F p = 45° (above) and the 30 GHz feed 
rotated by ipFP = 135° (below), both at position 1. In this figure the far field reference frame is rotated by (b) 
(j> XPD = -44.5° and (d) (j> X PD = -134.5° (see tabled). 



